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ABSTRACT 

Several studies in the recent past have inferred the existence of thermal inversions in some transit- 
ing hot Jupiter atmospheres. Given the limited data available, the inference of a thermal inversion 
depends critically on the chemical composition assumed for the atmosphere. In this study, we explore 
the degeneracies between thermal inversions and molecular abundances in four highly irradiated hot 
Jupiter atmospheres, day-side observations of which were previously reported to be consistent with 
thermal inversions based on Spitzer photometry. The four systems are: HD 209458b, HAT-P-7b, 
TrES-4, and TrES-2. We model the cxoplanet atmospheres using a 1-D line-by-line radiative transfer 
code with parametrized abundances and temperature structure, and with constraints of energy bal- 
ance and hydrostatic equilibrium. For each system, we explore the model parameter space with ~ 10 6 
models using a Markov chain Monte Carlo routine. Our results primarily suggest that a thorough 
exploration of the model parameter space is necessary to identify thermal inversions in hot Jupiter 
atmospheres. We find that existing observations of TrES-4 and TrES-2 can both be fit very precisely 
with models with and without thermal inversions, and with a wide range in chemical composition. On 
the other hand, observations of HD 209458b and HAT-P-7b are better fit with thermal inversions than 
without, as has been reported previously. Physically plausible non-inversion models of HD 209458b 
and HAT-P-7b fit the data only at the 1.7<7 observational errors; better fits require substantial en- 
hancement of methane and depletion of CO, which seems implausible in the very hot atmospheres 
considered here. Secondly, in the sample under consideration here, we do not see a correlation between 
irradiation levels and thermal inversions, given current data. Before JWST becomes available, near- 
IR observations from ground and with HST, along with existing Spitzer observations, can potentially 
resolve thermal inversions in some systems. Observations with only two channels of Warm Spitzer 
photometry and good S/N can likely identify or rule out thermal inversions if the difference between 
the fluxes in the 3.6 and 4.5 /jm channels is very high. 

Subject headings: planetary systems — planets and satellites: general — radiative transfer 



1. INTRODUCTION 

Observations of day-side atmospheres of several hot 
Jupiters have indicated the existence of thermal inver- 
sions. The inference of thermal inversions are a result 
of high S/N Spitzer observations (Knutson et al. 2008; 
Knutson et al. 2009; Machalek et al. 2009; O'Donovan 
et al. 2010; Christiansen et al. 2010) and concomitant 
theoretical modeling (Burrows et al. 2007, 2008; Fortney 
et al. 2008; Madhusudhan & Seager, 2009). Thermal in- 
versions have been reported based on flux excesses in cer- 
tain Spitzer channels. A natural explanation of channel- 
specific high flux is to invoke molecular emission features, 
as opposed to absorption features (Burrows et al. 2007, 
2008; Fortney et al. 2008). Under the assumption of lo- 
cal thermodynamic equilibrium (LTE), molecular emis- 
sion features form only in the presence of a thermal in- 
version, a region in the atmosphere where temperature 
increases outwards, as in the Earth's stratosphere. 

Preceding recent observations, early theoretical work 
on hot Jupiter atmospheres using self-consistent ID at- 
mosphere models predicted the existence of thermal 
inversions based on absorption due to TiO and VO 
(Hubeny et al. 2003, Fortney et al. 2006). More re- 
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cently, Fortney et al. 2008, classified hot Jupiters in two 
categories based on the degree of incident irradiation. 
The class of atmospheres with very high incident flux, 
dubbed "very hot Jupiters" or "pM" class, were consid- 
ered favorable to host thermal inversions caused due to 
gaseous TiO and VO absorption at low pressures, and at- 
mospheres with lower fluxes were predicted to be unlikely 
to host inversions owing to condensation of TiO/VO. The 
flux boundary for this dichotomy was, somewhat arbi- 
trarily, chosen to be ~ 10 9 ergs/s/cm 2 . However, while 
inferences of some recent observations have purported to 
violate this hypothesis (e.g. Machalek et al. 2009), oth- 
ers have found present observations insufficient to test 
this hypothesis (O'Donovan et al. 2010; Fressin et al. 
2010). 

The theory behind the absorbers causing thermal in- 
versions in hot Jupiters atmospheres is still under de- 
bate; discussed in detail in 12.11 Recent theoretical stud- 
ies suggest that TiO/VO may not be able to totally ac- 
count for the inferred thermal inversions (Spiegel et al. 
2009). Other alternatives proposed in recent studies in- 
clude strong UV/visible absorption due of photochemi- 
cally produced sulfur compounds (Zahnlc et al. 2009), 
and correlation of thermal inversions with chromosphcric 
activity of host stars (Knutson et al. 2010). Since the 
nature of absorbers causing the inversions is not known, 
models that have been successful in inferring thermal in- 
versions have either adopted a parametric absorber (Bur- 
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rows et al. 2008), or a parametrized temperature profile 
(see e.g. Madhusudhan & Seager, 2009). 

Existing models inferring thermal inversions rest on 
several assumptions and paramctrization. Models typ- 
ically invoke several free parameters to induce a ther- 
mal inversion in the temperature profile, as suggested by 
the data. The free parameters range from the location 
and magnitude of an unknown opacity source (e.g. Bur- 
rows et al. 2007 & 2008) to parametrizing the temper- 
ature profile itself (e.g. Madhusudhan & Seager, 2009). 
In addition, traditionally exoplanet atmosphere models 
span only a limited range in chemical abundances, of- 
ten assuming thcrmochcmical equilibrium (Barman et al. 
2005; Seager et al. 2005; Fortney et al. 2006; Burrows 
et al. 2007 & 2008). While some reported models re- 
quire thermal inversions to explain the observations, it is 
not known if the model parameter space has been thor- 
oughly exhausted. It remains to be seen whether some of 
the observations can be explained without thermal inver- 
sions, if we were to relax many of the assumptions and 
thoroughly explore the parameter space. 

In the present work, our primary goal is to understand 
the extent to which thermal inversions can be robustly 
inferred in hot Jupiter atmospheres, with current obser- 
vations. We choose a test sample of four systems for 
which Spitzer observations in the past have been reported 
to be consistent with thermal inversions. We then pur- 
sue a detailed exploration of the model parameter space 
to see the extent to which the observations can be ex- 
plained by models with and without thermal inversions. 
We accomplish this by computing large ensembles of in- 
version and non-inversion models (N ~ 10 6 ), explo ring 
the parameter space for best-fitting solutions. For each 
system, we report quantitatively how well the data can 
be fit by models with and without thermal inversions and 
with what ranges in atmospheric chemical composition. 
This approach also reveals the underlying correlations 
between the different chemical species and between the 
composition and temperature structure. And, since the 
systems considered here have different levels of irradia- 
tion, we also seek to understand if the presence or ab- 
sence of thermal inversions is correlated with the degree 
of irradiation, at the level of current data. 

We focus on four hot Jupiters for which existing ob- 
servations were reported to be consistent with thermal 
inversions in their day-side atmospheres, and for which 
photometric observations of thermal emission are avail- 
able in four or more Spitzer channels. The planets are: 
HD 209458b, HAT-P-7b, TrES-4, and TrES-2. Bur- 
rows et al. (2008) and Knutson et al. (2008) first re- 
ported a thermal inversion in HD 209458b, based on 
Spitzer photometry in five channels reported by Knut- 
son et al. (2008) and Deming et al. (2005). Madhusud- 
han & Seager (2009) confirmed the thermal inversion in 
HD 209458b, for model fits within the 1.4cr observational 
uncertainties. Additional observations of HD 209458b 
are available in the near-IR, obtained with HST NIC- 
MOS, but were explained by models with and without 
inversions alike (Swain et al. 2009). Knutson et al. 
(2009) reported a thermal inversion in TrES-4b based on 
five-channel Spitzer photometry using models based on 
Burrows et al. (2008). O 'Donovan et al. (2010) reported 
observations of TrES-2 in four Spitzer IRAC Channels, 
which were explained by models with and without ther- 



mal inversions; they preferred the inversion model which 
seemed more favorable amongst the set of models ex- 
plored in that study. Croll et al. (2010) reported a 
ground based detection of TrES-2 in the Ks band, and 
noted that their observation along with the four Spitzer 
observations could be explained equally well by mod- 
els with and without thermal inversions. More recently, 
however, Spiegel & Burrows (2010) reported that mod- 
els assuming radiative and chemical equilibrium cannot 
explain the Spitzer observations of TrES-2 without invok- 
ing thermal inversions. And, finally, Christiansen et al. 
(2010) reported a thermal inversion in HAT-P-7b, based 
on four Spitzer IRAC channels, along with an observa- 
tion in the visible by the Kepler Space Telescope (Borucki 
et al. 2009). While the Kepler point was not necessar- 
ily decisive in constraining the thermal inversion itself, 
it allowed constraints on the albedo, day-night redistri- 
bution and the TiO/VO abundance. Spiegel & Burrows 
(2010) confirmed the presence of a thermal inversion in 
HAT-P-7b using the observations of Christiansen et al. 
(2010). 

The observations of TrES-2 and HAT-P-7b noted 
above were first interpreted using the exoplanet at- 
mosphere model developed in Madhusudhan & Seager 
(2009), which is also used in the present work. In both 
those studies, i.e of O'Donovan et al. (2010) and Chris- 
tiansen et al. (2010), we had reported a representative 
set of models which explained the data. In the present 
study, we report a more exhaustive exploration of the 
model parameter space for these systems using a new 
parameter space exploration scheme described in Sec- 
tion o 

We discuss the theoretical and observational basis 
of thermal inversions in § 2. We explain the model 
framework in § 3, along with the parameter exploration 
method and the selection of systems for our study. In 
§ 4, we present our results, followed by a summary and 
discussion in § 5. 

2. THERMAL INVERSIONS IN HOT JUPITER 
ATMOSPHERES 

Several arguments have been proposed in the literature 
justifying the existence of thermal inversions in some hot 
Jupiter atmospheres. Compelling Spitzer IRAC observa- 
tions of some hot Jupiters suggest that anomalies in the 
form of planet-star flux excess in some channels cannot 
be explained without invoking thermal inversions (Bur- 
rows et al. 2007; Knutson et al. 2008; Madhusudhan 
& Seager, 2009; Swain et al. 2009). However, it is not 
known whether such an inference is an outcome of the 
model input choices used to infer the observations. On 
the other hand, independent of observations, several the- 
oretical arguments support the existence of thermal in- 
versions in some hot Jupiter atmospheres (Hubeny et al. 
2003; Burrows et al. 2007; Fortney et al. 2008). 

In this section, we explore the arguments in favor of 
thermal inversions, and motivate the framework under 
which they can be tested. We begin with the theoretical 
motivation for thermal inversions, followed by arguments 
leading to inference of thermal inversions from observa- 
tions. We then pose the question of whether the obser- 
vations can be explained without thermal inversions if 
some of the model constraints are relaxed. Finally, we 
set up the framework in which the requirement of ther- 
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mal inversions can be robustly tested. 

2.1. Theoretical Basis for Thermal Inversions 

Thermal inversions are a natural consequence of visi- 
ble/UV absorption of incident star light high in the atmo- 
sphere. For an isolated planetary atmosphere in hydro- 
static equilibrium, and no local energy sources, the atmo- 
spheric temperature decreases with pressure (i.e with in- 
creasing distance from center) ; the atmosphere is heated 
from below, and cools monotonically outwards. How- 
ever, in planetary atmospheres irradiated by the host 
star, strong optical/UV absorbers in the higher layers of 
the atmosphere can intercept part of the incident star 
light. Such local deposition of energy results in a zone in 
the planetary atmosphere where temperature increases 
outward, i.e a "thermal inversion" . Most solar system 
planets have thermal inversions in their atmospheres. In 
Earth's atmosphere, for example, a thermal inversion is 
caused by ozone (O3), which is a strong absorber in the 
UV (Chamberlain, 1978). And, in Jupiter's atmosphere, 
a thermal inversion is caused by strong absorption in the 
visible by haze resulting from methane photochemistry. 
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Fig. 1. — Illustration of P-T profiles. The red and blue curves 
show P-T profiles with and without a thermal inversion, respec- 
tively. The molecular features due to each of these profiles are 
shown in Figure [2] 

The thermal inversion causing absorbers of solar sys- 
tem planet atmospheres do not survive the temperatures 
of hot Jupiters. Nevertheless, it has been proposed that 



thermal inversions could be formed in the atmospheres 
of very hot Jupiters due to strong absorption of incident 
stellar radiation in the visible by gaseous TiO and VO 
(Hubeny et al. 2003; Burrows et al. 2007; Fortney et 
al. 2008). However, Spiegel et al. (2009) suggested that 
at mbar pressures where thermal inversions are required 
to explain the observations, TiO and VO may not be 
abundant in the required amounts. TiO, being a heavy 
species, requires a substantial amount of vertical mixing 
to keep it aloft; the K zz required is 10 7 — 10 11 cm 2 /s, 
for 0.1-10 /im condensate size, if a cold trap is present. 
Since neither the particle sizes nor K zz are known, it is 
uncertain if TiO might be able to explain thermal inver- 
sions. And, VO is unlikely to help owing to its lower (so- 
lar) abundance and lower visible absorption compared to 
TiO (Spiegel et al. 2009). Zahnle et al. (2009) reported 
that photochcmically produced Sulfur compounds, HS 
and S2, could have high UV and visible opacities at the 
temperatures relevant to hot jupiter atmospheres, mak- 
ing them potential candidates for causing thermal in- 
versions in hot jupiters. More recently, Knutson et al. 
(2010) suggested that the presence of a thermal inver- 
sion could be inversely correlated with the activity level 
of the host star, UV flux from the more active stars po- 
tentially destroying inversion-causing absorbers. Despite 
the continuing debate on the inversion-causing absorbers 
in hot jupiter atmospheres, thermal inversions have been 
invoked typically by either adopting a parametric opac- 
ity source of unknown nature (Burrows et al. 2008), or 
by parametrizing the temperature profile (see e.g. Mad- 
husudhan & Seager, 2009). 

2.2. Review of Observational Inference of Thermal 
Inversions 

The inferences of thermal inversions are motivated by 
observations of dayside atmospheres of transiting hot 
Jupiters using Spitzer photometry. Observations of some 
hot Jupiters indicate excess emission in some Spitzer 
channels over others. For example, in the four IRAC 
observations of HD 209458b reported by Knutson et 
al.(2008), there is indication of excess emission in the 
4.5 /im and 5.8 /im channels. The observations show 
a markedly higher planet-star flux ratio in the 4.5 /im 
channel compared to the neighboring 3.6 /im channel. 
And, the flux ratio in the 5.8 /im channel is higher than 
that in the adjacent channels at 4.5 /im and 8 /im. 

The thermal emission spectrum of a planet is influ- 
enced by a combination of the atmospheric temperature 
structure and molecular absorption. Let us consider a 
hypothetical hot Jupiter atmosphere with the stellar and 
planetary properties of HD 209458b, and consisting of 
gaseous H2O, CO, CH4, and CO2; at nominal mixing ra- 
tios, close to assumption of chemical equilibrium at solar 
abundances. In this particular case, we assume that the 
molecular species are all well mixed in the atmosphere, 
i.e. uniform volume mixing ratios over the whole atmo- 
sphere. Figure[2]shows the spectral features of each of the 
molecules in such an atmosphere, for P-T profiles with 
and without a thermal inversion (shown in Figure [1}. 
The assumed mixing ratios of the molecules are: H2O = 
10" 4 , CO = 10- 3 , CH 4 = 10- 5 , C0 2 = 5x 10~ 7 . As is 
demonstrated in Figure [U the atmosphere with a ther- 
mal inversion gives rise to molecular emission features, 
whereas the one without a thermal inversion has absorp- 
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Fig. 2. — Illustration of molecular features. Each panel shows molecular line features due to a single molecule in a hypothetical atmosphere. 
The stellar and planetary properties are assumed to be those of the HD 209458b system. The blue curves show the spectral features of an 
atmosphere without a thermal inversion, and red curves show features of an atmosphere with a thermal inversion. The corresponding P-T 
profiles arc shown in Figurc[T| The mixing ratios, with respect to H2, in the four panels are assumed to be: H2O = 10 -4 , CO = 10 — 3 , 
CH4 = 10 — 5 , CO2 = 5 X 10 . The continuum in each spectrum is of the black-body of the bottom of the atmosphere with features due 
to H2-H2 collision induced opacities. The green curves show the six Spitzer band passes at 3.6, 4.5, 5.8, 8, 16, and 24 (mi. 



tion features. 

The excess emission in some Spitzer channels over oth- 
ers can be qualitatively explained by a thermal inversion 
along with some key molecular features. As can be seen 
from Figure [2l H 2 has several spectral features in the 
3.6 /im, 5.8 /Urn, 8 fim, 16 /im, and 24 fan Spitzer chan- 
nels. CH4 has strong features almost exclusively in the 
3.6 /xm and 8 /im channels. CO has a strong feature in 
the 4.5 /im channel, also contributing to the 5.8 /im chan- 
nel. And, CO2 has strong features in the 4.5 /im and 16 
/im channels. The 4.5 /im feature of CO2 is degenerate 
with the contribution of CO in the same channel. The 
high fluxes in the 4.5 /im and 5.8 /im channels, as seen 
in HD 209458b for example, can therefore be explained 
simply by having strong emission due to CO and CO2, 
and only moderate emission due to CH4 and H2O (Mad- 
husudhan & Seager, 2009). And, since emission features 
can form only due to a thermal inversion, the observa- 
tions can be interpreted as suggesting the presence of a 
thermal inversion along with CO and/or C02- 

The molecular species required to explain the obser- 
vations arc physically plausible. In chemical equilibrium 
(Burrows & Sharp, 1999), CO occurs naturally in the 
temperature range of HD 209458b (see P-T profiles in 
Burrows et al. 2008, Madhusudhan & Seager, 2009, 



Swain et al. 2009). And, CO2 mixing ratios up to ~ 10 6 
arc allowed by equilibrium chemistry and/or photochem- 
istry (Liang et al. 2003; Zahnlc et al. 2009). More 
generally, hot Jupiters for which thermal inversions are 
predicted, are characterized by their very hot day-side 
atmospheres (Fortney et al. 2008), suggesting that dom- 
inant contribution to the emergent spectra is expected 
from CO (Burrows & Sharp, 1999). 

Care must be exercised while exploring the space of at- 
mospheric composition in order to fit the data. Since a 
high flux in the 4.5 /im channel can be explained by CO 
and/or CO2, it is possible for a fitting model to infer the 
lack of CO in the atmosphere, by allowing an implausibly 
high C02- Such a proposition would likely be implausible 
given the hot temperatures where CO is expected to be 
the dominant form of carbon. The degeneracy between 
CO and CO2 can be broken by fitting the model to the 
16 /im Spitzer IRS photometry, where available. The 16 
/im channel has dominant contribution due to a strong 
feature of CO2 alone. For example, an observation of 
HD 209458b in the 16 /im channel by Deming (personal 
communication, 2009) was used in Madhusudhan & Sea- 
ger (2009) to place simultaneous constraints on CO and 
C02- The same observation is also used in the present 
work. 
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2.3. Alternate Qualitative Interpretation 

Is there any conceivable scenario, independent of ex- 
isting models in which the observations can be explained 
without thermal inversions? A few qualitative alterna- 
tives seem feasible. Let us re-consider the situation ex- 
plained in § 12.21 where the fluxes in the 4.5 /im and 5.8 
jitm IRAC channels arc higher than those in the 3.6 /an 
and 8 /im channels, respectively. Let us now investigate 
if we can explain the same observations with a planet at- 
mosphere which has no thermal inversion (the blue model 
in Figure Q] and Figure [5]). In this case, instead of con- 
sidering the emission features in the 4.5 //m and 5.8 /jm 
channels, one can consider the absorption features in the 
3.6 fim and 8 /jm channels, to explain the same flux dif- 
ferential between adjacent channels. The fitting model 
would then require strong absorption features in the 3.6 
/zm and 8 /zm channels, and weaker absorption in the 
4.5 /jm and 5.8 zim channels. One conceivable solution is 
provided by methane (CH4) which has strong absorption 
features only in the 3.6 /im and 8 /im channels. So, in 
principle, a high contribution due to CH4, over CO or 
CO2, could provide the required absorption signatures. 
Although H2O also has features in the 8 (im channel, a 
high abundance of it may not be desirable since it would 
lead to absorption in the 5.8 /im channel. In presenting 
this solution, we have exploited the degeneracy between 
the presence or absence of a thermal inversion and com- 
plementary molecular features. 

Despite our simple qualitative explanation of a 
non-thcrmal-inversion fit to observations like those of 
HD 209458b, whether or not a model without inversion 
fits the data, and is physically plausible, is subject to 
further tests. Firstly, CH4 has stronger absorption in 
the 8 /jm channel than in the 3.6 /jm channel. In other 
words, for the same molecular composition of CH4, the 
absorption in the 8 /jm channel with respect to the 5.8 
/im channel can be deeper than the absorption in the 3.6 
/xm channel with respect to the 4.5 /zm channel. This is 
contrary to what is required by the IRAC observations 
of HD 209458b (Knutson et al. 2008). 

A fitting model without a thermal inversion must also 
satisfy the constraint of energy balance. Explaining the 
high flux in the 4.5 /jm channel with a non-thermal in- 
version model means that the black-body continuum of 
the spectrum must be at the level of the 4.5 /im point or 
higher. While such a high emergent flux balances the in- 
cident stellar flux remains to be verified. Therefore, it is 
not certain that a non-inversion configuration which fits 
the data necessarily satisfies the fundamental constraint 
of energy balance. This latter point could, in principle, 
be obviated by P-T profiles steep enough to produce deep 
spectral features in most parts of the spectrum except in 
the 4.5/im and 5.8 /im channels. Finally, even if the 
model fits and maintains energy balance, it is not clear 
how such extremely hot atmospheres can be dominated 
by CH4 over CO. Future theoretical work might explain 
this possibility. 

2.4. A Test for Thermal Inversions 

A rigorous inference of a thermal inversions from a 
given set of observations would involve running a large 
population of models thoroughly exploring the parame- 
ter space in search of non-inversion models that fit the 



data. Not finding a statistically significant fit in such an 
exploration would constitute strong evidence in favor of 
thermal inversions from the data set in question. In a 
recent work, we demonstrated the capability of running 
~ 10 7 1-D models on a parameter grid with manage- 
able computational resources (Madhusudhan & Seager, 
2009). In the same study, we had explored models on 
predetermined grids in the parameter space of models 
for HD 189733b and HD 209458b. We were then able 
to constrain the extent of possible thermal inversions in 
HD 209458b that fit the Spitzer broadband photometry. 
Such a capability was possible because of the efficient 
paramctrization of the model temperature structure and 
molecular abundances. 

In this study, we combine the model developed in Mad- 
husudhan & Seager (2009) with an efficient parameter 
space exploration procedure to test the requirement of 
thermal inversions for a select sample of hot Jupiters. We 
consider four hot Jupiters at different levels of irradiation 
for which observations have been known to be consistent 
with thermal inversions in their respective atmospheres, 
and for which Spitzer photometry is available in four or 
more channels. We run ~ 10 6 models, with and without 
thermal inversions, for each planet under consideration, 
and report goodness-of-fit contours in the space of atmo- 
spheric composition and temperature structure. In what 
follows, we describe our model set-up, the optimization 
algorithm, and the systems considered in this study. 

3. MODEL AND METHOD 
3.1. Model 

Our model atmosphere includes a ID parametric P-T 
profile coupled with line-by-line radiative transfer, hydro- 
static equilibrium, and the requirement of energy balance 
at the top of the atmosphere (Madhusudhan & Seager, 
2009). We consider 100 atmospheric layers in the pres- 
sure range between 10~ 5 — 100 bar. The key aspect of 
our model is the paramctrization of the P-T profile and 
the chemical composition, which allows us to run large 
ensembles of models, exploring the parameter space, in 
a computationally efficient manner. 

The major difference of our model from traditional at- 
mosphere models is in the treatment of energy balance. 
Our model requires energy balance at the top of the at- 
mosphere, instead of an iterative scheme to ensure layer- 
by-layer radiative (or radiative + convective) equilibrium 
as is done in conventional models. We note that the 
requirement of layer- by- layer radiative equilibrium in a 
1-D model is not strictly physical since complex hydrody- 
namics flows in highly irradiated hot Jupiter atmospheres 
can alter the temperature structure away from radiative 
equilibrium (Showman et al. 2009). The global energy 
balance, e.g., at the top of the atmosphere, however, is a 
strict requirement. For a given set of model parameters, 
we require that the net energy output at the top of the 
atmosphere is less than or equal to the net energy input 
due to the incident stellar flux; a deficit indicates energy 
redistributed to the night-side. Models where the emer- 
gent flux is greater than the incident flux are discarded 
(see Madhusudhan & Seager, 2009). By running a large 
number of (~ 10 7 ) models in the parameter space, and 
discarding those that did not satisfy energy balance, we 
were left with a population of models that satisfied en- 
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ergy balance. 

We parameterize the chemical abundances of the 
molecular species by considering deviations over chem- 
ical equilibrium (explained at length in Madhusudhan 
and Seager, 2009). For each molecule under considera- 
tion, we compute its mixing ratio in a layer of the atmo- 
sphere by multiplying a parametric factor to the mixing 
ratio that would be expected under thermochemical equi- 
librium with solar abundances (TEq). The parametric 
factor for a given molecule is constant over the entire at- 
mosphere, i.e the mixing ratio profile of the molecule over 
the entire atmosphere would be shifted relative to that 
obtained from TE Q by the constant factor. We use this 
treatment for H2O, CO, and CH4. For CO2, we perturb 
over a uniform mixing ratio of 10~ 6 (1 ppmv or 10~ 6 is 
just a reference; it is the approximate mixing ratio of CO2 
expected from TEq and/or photochemistry, for 5 x solar 
metallicity and T ~ 2000A"; Zahnlc ct al. 2009, Liang 
et al. 2003). Thus corresponding to the four prominent 
molecules, we have four parameters: /h 2 Oj /co, /ch 4 , 
and /co 2 • We reiterate that fx is not the absolute mix- 
ing ratio of "X" . It is the ratio between the mixing ratio 
of "X" and the mixing ratio under TE Q ; except for CO2, 
for which it is with respect to 10~ 6 . Our models also 
include NH 3 , fixed at the TE value, and TiO and VO 
at solar abundances. Additionally, we include H 2 -H 2 col- 
lison induced cross-sections, which are a source of con- 
tinuum opacity. Our H2O, CH4, CO and NH3 molecular 
line data are from Freedman et al. (2008), and refer- 
ences therein. Our CO2 data are from Freedman (per- 
sonal communication) and Rothman et al. (2005). And, 
we obtain the H2-H2 collision-induced opacities from Bo- 
rysow et al. (1997), and Borysow (2002). We use a Ku- 
rucz model for the stellar spectrum (Castelli & Kurucz, 
2004). 

In the current work, we have made one significant 
change to our approach in Madhusudhan & Seager 
(2009). Previously we had run models on a predeter- 
mined grid, chosen based on some model independent 
arguments outlined in that work. While we were able 
to run a large population (tens of millions) of models in 
that approach, the grid resolution was still coarse and 
sampled evenly over evidently unnecessary regions of the 
parameter space. In this approach, we run the models 
using a more efficient parameter space exploration pro- 
cedure, allowing us to sample the desired error surfaces 
at much higher resolution. 

3.2. Parameter Space Exploration 

Our primary requirement in this work is to be able 
to explore the model parameter space at fine resolution. 
Even a single plausible model without thermal inversions, 
in a million models, would still be evidence against the 
requirement of thermal inversion by a given set of obser- 
vations. 

We use the Markov chain Monte Carlo (MCMC) 
method to explore the parameter space of models without 
thermal inversions. The MCMC method is a Bayesian 
parameter estimation algorithm which allows the calcu- 
lation of posterior probability distributions of the model 
parameters conditional to a given set of observations. An 
extensive body of literature exists on the applications 
of MCMC for parameter estimation (Gilks ct al. 1998; 
Tegmark ct al. 2004; Ford et al. 2005). The MCMC 



method allows an efficient means of exploring the param- 
eter space in search of a global solution, with very fine 
sampling in the allowed range of parameter values. In 
this work, however, the observations are always less than 
the number of parameters, i.e there is no unique solu- 
tion. However, it is still possible to explore the parameter 
space and find contours in the error surface of some mea- 
sure of fit. We, therefore, use the MCMC method with 
a Metropolis-Hastings scheme within the Gibbs sampler, 
for fine sampling of the model parameter space. 

Our model described in § l3.1l above has ten free param- 
eters. Six parameters concern the P-T profile: To, Pi, 
P 2 , P3, ai, and «2 (Madhusudhan & Seager, 2009). And, 
four parameters correspond to the departures of molecu- 
lar abundances from the reference abundances described 
in § EIU /h 2 o, /co, /ch 4 , and / C o 2 - 

We define some physically motivated boundaries in the 
parameter space explored by the Markov chain. We 
impose the constraint of global energy balance by re- 
stricting r] to [0.0,1.0], where, 77 = (1 — A)(l — f r ) is 
the ratio of emergent flux output on the day-side to in- 
cident stellar flux input on the day-side, weighted ap- 
propriately (Madhusudhan & Seager, 2009). Here, A 
is the Bond Albedo and f r is the day-night energy re- 
distribution. And, we impose some nominal bound- 
aries on the temperatures and departures from equilib- 
rium chemistry. We explore a wide range of deviations 
from chemical equilibrium (Burrows & Sharp, 1999), 
empirically selected so as to be general enough. For 
models without thermal inversions, we set the bound- 
aries as -9 < log(/ H2 o) < 9, -9 < log(/co) < 5, 
-5 < log(/ C H 4 ) < 10, -5 < log(/ c02 ) < 4. The lim- 
its are similar for models with thermal inversions, except 
the lower boundaries for log(/H 2 o) and log(/co), which 
are set at -4 and -5, respectively. We report all those 
models which have the overall elemental C/H and O/H 
abundances within the broad range of (10~ 2 — I0 2 ) x 
solar (Anders & Grevesse, 1989; Burrows & Sharp, 1999; 
but c.f. Asplund & Grevesse, 2005; Allende-Prieto et al. 
2002). For the temperature structure, the constraint 
of no thermal inversion is imposed by requiring that 
Pi > Pz- The "fit" parameters for the MCMC arc T , 
log(Pi), log(P 2 ), log(P 3 ), ai, a 2 , k>g(/ H2 o), log(/co), 
Iog(/cH 4 ); and log(/co 2 )- We consider uniform priors in 
all the parameters. For each system under consideration, 
we run one chain of 10 6 links for models with thermal in- 
version and one for models without thermal inversion. 
Our parametric P — T profile provides a simple means to 
demarcate between inversion and non-inversion models. 
The condition for the P-T profile to have no thermal in- 
version is P2 < P\ < P3. And, that to have a thermal 
inversion is Pi < P2 < P3. 

3.3. Quantitative Measure of Fit 

Central to our analysis is the definition of what con- 
stitutes a "fit" to the data. We can only report to what 
extent a model fits the data relative to the "observa- 
tional" uncertainties, i.e within the 1 a error bars, or 1.5 
a errors, and so on. Given that the number of broadband 
observations are typically less than the number of model 
parameters, we cannot report a formal fit with confidence 
levels. Nevertheless, we evaluate our models based on the 
£ 2 statistic, defined as x 2 /Nobs (Madhusudhan & Seager, 
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Fig. 3. — Pressure- Temperatures (P-T) profiles for the four systems. Each panel shows the P-T profiles with no thermal inversions 
that explain the observations at different levels of fit. Best-fitting profiles with thermal inversions for each of these systems have been 
reported in the literature referred in the text. For HD 209458b. the profiles in magenta correspond to models that fit the observations with 
2 < § 2 < 2.25; the best-fit non-inversion model has £ 2 = 2.04. For HAT-P-7b, the red profiles correspond to models fitting the observations 
to within 1.6 < £ 2 < 2; the best-fit model had £ 2 = 1.65. The brown profiles for both HD 209458b and HAT-P-7b are 30 profiles that 
fit best, shown for illustration. For TrES-4 and TrES-2 the red profiles correspond to models that fit within 1 < £ 2 < 2, and the purple 
profiles fit to within f; 2 < 1; only 100 randomly chosen profiles from each category are shown, for clarity. 
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Fig. 4. — Departures from chemical equili briu m of CH4 and CO. The dots indicate the regions in the space of CO and CH4 mixing ratios 
explored by the MCMC chain (see Section l3.2| h each dot is a model realization. All m odels with C/H and O/H within (10 — 2 — 10 2 ) X 
solar are shown. The boundaries in the composition space are described in section [372] For each planet, f^o and fcH 4 are the departures 
in the mixing ratios of CO and C H4 fr om those corresponding to thermochemical equilibrium with solar abundances (TEq) for the same 
temperature structure (see Scction l3.1j l, For example, fco = 1 implies a CO concentration that is at TEq. The left (right) panel for each 
system shows constraints on models without (with) thermal inversions. Non-inversion and inversion models are labelled with (Nl) and (I), 
respectively. The purple, orange, green and black colors correspond to J 2 < 1, 1 < £ 2 < 2, 2 < § 2 < 3, and 3 < £ 2 < 4, respectively. The 
blue dots for non-inversion models in HD 209458b and HAT-P-7b correspond to 2.0 < £ 2 < 2.25; the best fitting model for HD 209458b 
has a £ 2 = 2.04, and that for HAT-P-7b has £ 2 = 1.65. 
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2009): 

^2 ^ ^ y f fi,model fi.obs j 

Nobs -, \ &i.obs J 

i— 1 v 7 

where, fi, m odei and /i, b s are the model and observed flux 
ratios, respectively, and <7i. b s is the 1 c measurement 
uncertainty. N D bs is the number of observations. 

For each system, we report the best value of £ 2 we 
find with a non-inversion model. We also present the 
range in parameter space which fit the observations at 
different levels of £ 2 , for example, £ 2 < 1, £ 2 < 2, or 
higher, as applicable. In this framework, £ 2 < 1 means 
the model fits the observations to within the 1 a error 
bars on average. And, a £ 2 < 2 means a fit to within 
\/2 = 1.4 a of the error bars on average, and so on. 
4. RESULTS 

In this section, we report the constraints on the chemi- 
cal compositions of the day-side atmospheres of the plan- 
ets in our study. For each planet, we present the range 
of composition, temperature structure, and day-night en- 
ergy redistribution required by the best-fit models with 
no thermal inversions. We also present the constraints 
on the composition and day-night energy redistribution 
of models with thermal inversions for each system. 

4.1. HD 209458b 

We consider planet-star flux contrasts of HD 209458b 
in six channels of Spitzer broadband photometry. The 
data include four IRAC observations reported by Knut- 
son et al. (2008), and observations in the 16 /xm IRS 
channel and the 24 /im MIPS channel, by Deming (per- 
sonal communication) and Deming et al.(2005), respec- 
tively. We focus on these observations which were re- 
ported in the literature as suggestive of a thermal inver- 
sion in HD 209458b, especially the four IRAC observa- 
tions. 

Our results indicate that HD 209458b is a likely can- 
didate to host a thermal inversion in its day-side atmo- 
sphere. However, whether or not HD 209458b actually 
has a thermal inversion depends on the level of fit, and 
the physical plausibility of the fitting models one is will- 
ing to consider. Figure [3] shows populations of pressure- 
temperature (P-T) profiles with no thermal inversions 
which fit the observations at different levels of £ 2 . The 
corresponding constraints on the atmospheric composi- 
tion are shown in Table [IJ and the constraints for inver- 
sion models are shown in Tabled Figure @] and Figure [S] 
show departures of molecular species from TEq , for mod- 
els of HD 209458b with and without thermal inversions. 

The observations require a thermal inversion in the at- 
mosphere of HD 209458b at the £ 2 = 2 level. The best 
fitting model with no thermal inversion has a £ 2 of 2.04, 
implying a fit at 1.43 a (i.e V2.04) of the observations, on 
average. And, even at this level of fit, the models show 
substantial departures from thermochemical equilibrium 
assuming solar abundances (TE ). Figure 0] shows the 
departures in the mixing ratios of CO and CH4 from TEq 
at different levels of fit. At the 2.0 < £ 2 < 2.25 surface 
(shown in blue dots for the non-inversion case) , it can be 
seen that non-inversion models of HD 209458b require a 
depletion of CO of at least 10 -2 times TEq (the depar- 
tures shown in Figure [4] are in fraction with respect to 



TEq). This low a mixing ratio of CO can, in principle, be 
achieved by having similarly low abundances of C and O 
relative to solar. However, the simultaneous requirement 
of an overabundance of CH4, is hard to explain. Thus, 
non-inversion models fitting the observations at the 1.5 
a errors seem physically implausible. 

The observations can be explained by physically plau- 
sible non-inversion models at the £ 2 ~ 3 level, meaning 
a fit at the 1.7 a observational errors. As shown in Fig- 
ure 21 the £ 2 = 3 level (region with green dots) allows for 
non-inversion models which have CO within a factor of 
~ 10 from TEq. Such small a factor can potentially be 
explained either by just having different C/H and O/H 
abundances or due to non-equilibrium processes (Cooper 
& Showman, 2006; Zahnle et al. 2009; Line et al. 2010; 
Madhusudhan & Seager, 2010). 

The observations can be fit to within £ 2 = 1 by models 
with thermal inversions, as has been demonstrated previ- 
ously by Madhusudhan & Seager (2009). The molecular 
mixing ratios for inversion models as constrained by the 
observations are shown in Table [2] Figures |4] and [5] show 
the deviations from TEq required by the inversion mod- 
els. We find that the best-fitting models, within £ 2 = 1, 
allow compositions with CO and CH4 deviant from TEq 
by a factor of ~ 10 and higher. However, if we consider 
the £ 2 = 2 surface, the observations can be fit with in- 
version models containing close to TE values of CO and 
CH4, with C/H and O/H abundances only slightly en- 
hanced over solar. One potential problem in the inversion 
scenario, however, is the requirement of low H2O by the 
best-fit models. The low observed flux in the 3 jim and 
24 /im channels requires low H2O, which is contrary to 
the high H2O requirement imposed by the high observed 
flux in the 5.8 /im channel. This problem has been pre- 
viously discussed in the literature (Deming et al. 2005, 
Seager et al. 2005, Madhusudhan & Seager, 2009), and 
is a subject for future studies. 

Two best-fit model spectra, with and without a ther- 
mal inversion, for HD 209458b are shown in Figure El 
along with the P-T profiles. The corresponding at- 
mospheric composition and day-night redistribution are 
shown in Table [3] As can be seen from Figure |H the 
dominant source of error for the non-inversion model 
comes from the high flux in the 5.8 /im IRAC obser- 
vation. While most of the observations can be fit at the 
~ lcr error bars, the 5.8 /im is fit only at ~ 2.5er. The 
dominant source of opacity in this channel comes from 
H2O. In the non-inversion scenario, a low water content 
could potentially explain the lack of absorption in this 
channel but will over predict the flux in the 24 /im chan- 
nel. 

4.2. EAT-P-lb 

HAT-P-7b is one of the hottest transiting hot Jupiters 
known. Being the hottest of our sample of planets, it 
is also the most expected hot Jupiter in our sample to 
host a thermal inversion. We use the day-side observa- 
tions of HAT-P-7b reported in four channels of Spitzer 
IRAC photometry by Christiansen et al. (2010). The 
constraints on the molecular abundances, and departures 
from equilibrium, for the inversion and non-inversion 
models are shown in Table [TJ Table [31 Figure 0] and 
Figure [5l 

Our results confirm previous findings that the inver- 
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TABLE 1 

Constraints on the atmospheric properties for models without thermal inversions 





HD 209458b a 


HAT-P-7b a 




TrES-4 






TrES-2 






2 < £ 2 < 2.25 


1.5 < £ 2 < 2 


e 


< 1 


e 


< 2 


e 


< 1 


e 


< 2 


H 2 O b 


1CT 5 - 2 x 10" 3 


10~ 12 - 2 x 10- 5 


10- 12 - 


2 x 10" 3 


10- 13 


- 0.15 


8 x 10" 


6 - 0.16 


10- 12 


- 0.16 


CO 


lCT 12 - 10" 5 


3 x 10~ 8 - 10~ 4 


10- 13 


- 0.06 


10- 13 


- 0.07 


10- 13 


- 0.04 


10- 13 


- 0.07 


CH 4 


6 x 1CT 6 - 2 x 10" 4 


2 x 10~ 3 - 0.05 


10~ 15 


- 0.07 


10- 15 


- 0.07 


10" 15 


- 0.07 


10- 18 


- 0.07 


C0 2 


10- 11 - 4 x icr 8 


9 x 10~ 7 - 9 x 10~ 5 


3 x 10" 


7 - 0.02 


10- 12 


- 0.03 


10- 11 


- 10" 5 


10- 12 - 


7 x 10~ 3 


C/O 


0.03 - 0.82 


15 - 2 x 10 3 


0.44- 


4 x 10 3 


10- 4 - 


4 x 10 3 


5 x 10~ 5 


- 4 x 10 3 


4 x 10" 5 


- 4 x 10 3 


rf 


0.5 - 1.0 


0.7 - 1.0 


0.37 


- 1.0 


0.33 


- 1.0 


0.36 


- 1.0 


0.29 


- 1.0 



a For HD 209458b and HAT-P-7b, the best-fit non-inversion model has £ of 2.04 and 1.65, respectively, and hence the reported ranges of £ . 
Sec text for details. 

k The molecular mixing ratios are quoted as ratios by number with respect to molecular hydrogen. 

c rj — (1 — / r )(l — A), where, f r is the day-night redistribution, and A is the bond albedo. (l-?7) gives the maximum day-night redistribution 
allowed by the model, i.e assuming zero albedo. 



TABLE 2 

Constraints on the atmospheric properties for models with thermal inversions 







HD 209458b 






HAT-P-7b 




TrES-4 a 


TrES-2 a 


e 


< 1 


e 


< 2 


e 


< 1 


e 


< 2 


€ 2 < 1 


i 2 < 1 


H 2 O b 


5 x 10~ 8 


-3 x 10~ 6 


4 x 10" 


8 _ 10 -5 


10- 4 


- 0.16 


3 x 10" 


7 - 0.16 


5 x 10~ 8 - 0.13 


5 x 10~ 8 - 0.16 


CO 


3 x 10" 


4 - 0.06 


3 x 10" 


4 - 0.07 


3 x 10" 


9 -0.07 


IO- 9 


- 0.07 


4 x 10~ 9 - 0.07 


10~ 9 - 0.07 


CH 4 


9 x 10- 


4 - 0.06 


10- 11 


- 0.06 


10- 14 


- 0.07 


10- 15 


- 0.07 


10~ 15 - 0.07 


10~ 14 - 0.07 


C0 2 


Hr 11 - 


5 x 10~ 7 


10- 12 - 


8 x KT 7 


10- 11 - 


2 x 10~ 4 


10- 12 


- 0.01 


10~ n - 0.02 


10" 11 - 0.02 


C/O 


1 - 


■ 146 


1 - 


105 


5 x 10 


" 5 - 53 


5 x 10" 


" 8 - 207 


4 x 10~ 4 - 2 x 10 3 


6 x 10" 5 - 10 3 


rf 


0.27 


-0.60 


0.25 


- 0.84 


0.49 


- 0.83 


0.41 


- 0.85 


0.36 - 1.0 


0.23 - 1.0 



For TrES-2 and TrES-4, the constraints at the £ — 2 level are almost identical to those at the £ — 1 level, and hence we do not report 
them here. 

b The molecular mixing ratios are quoted as ratios by number with respect to molecular hydrogen. 

c 77 — (1 — / r )(l — A), where, f r is the day-night redistribution, and A is the bond albedo, (l-ij) gives the maximum day-night redistribution 
allowed by the model, i.e assuming zero albedo. 



sions models fit the observations of HAT-P-7b better 
than models without thermal inversions. As shown in 
Figure 21 inversion models can fit the observations to 
within the 1-cr errors (i.e. £ 2 < 1) for a wide range of 
methane and CO concentrations, including those close 
to TE© values. On the other hand, the best fitting non- 
inversion model has a £ 2 = 1.65, indicating a fit at 1.3-cr 
errors. Even then, the best fitting non-inversion mod- 
els of HAT-P-7b shown in Figure |4] (blue dots) require 
methane abundances that are over five orders of magni- 
tude greater than the TE© values, which are seemingly 
implausible. 

The observations of HAT-P-7b can plausibly be ex- 
plained without a thermal inversion at the fj 2 = 3 level, 
i.e fits at 1.7-er errors (green dots). At this level, mod- 
els with methane and CO concentrations only marginally 
deviant from TE© values can explain the data. However, 
such a degree of fit may not be statistically representative 
of the true nature of the planet atmosphere. A sample 
of non-inversion temperature profiles at different levels 
of fit is shown in Figure H The P-T profiles at £ 2 < 2 
are shown in Figure [3] The brown profiles show the 30 
best-fit P-T profiles below £ 2 = 2. 

Two model spectra, corresponding to models with and 
without a thermal inversion, are shown in FigureJSl along 
with the corresponding temperature profiles. The model 
parameters for each case are shown in Table [3] As shown 
in Figure |6l the non-inversion model is unable to fit the 



low flux at 3.6 /jm and the high flux at 8 /^m, simul- 
taneously. The reason for this behavior is because both 
channels have absorption features due to methane. A low 
(high) methane concentration in a non-inversion model 
causes high (low) observed flux in both the channels, con- 
trary to what is observed. On the other hand, the inver- 
sion model easily explains all the observations. The low 
flux in the 3.6 micron channel is explained by the tem- 
perature decreasing outward in the lower layers, and the 
high flux at 8 microns is explained by the high tempera- 
tures, due to the thermal inversion, in the higher layers 
where the contribution to the 8 /an channel peaks. 

The day-night redistribution is well-constrained by the 
data. At the £ 2 = 2 level, the constraint on 77 — 
(1 — A)(l — f r ) is 0.7 - 1.0, for non-inversion models. 
This range allows for a maximum redistribution of 0.26, 
assuming zero albedo, and hence implies relatively ineffi- 
cient advection of energy to the night side. This finding 
is consistent with the low redistribution of an inversion 
model reported in Christiansen et al. (2010), and the 
finding of low redistribution in the visible light curve of 
Kepler. The inversion models on the other hand allow 
for a wider range of r] = 0.41 — 0.85 at the £ 2 < 2 level, al- 
lowing for models with efficient day-night redistribution 
and those otherwise. 

In principle, the allowed parameter space of the models 
can be further restricted by the observation in the visi- 
ble (centered at 0.63 /_tm) obtained by the Kepler space 
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Fig. 6. — Sample model spectra for each system. Two model spectra, corresponding to models with and without a thermal inversion, are 
presented for each system. The models represent a balance between degree of fit and physical plausibility. The corresponding P-T profiles 
are shown in the insets. The atmospheric compositions and day-night redistribution corresponding to each model are shown in Table [3] In 
each panel, the black circles with error bars (and the upper-limit for TrES-4) are the available Spitzer observations (see text for details). 
The red and green circles are the channel integrated model spectra in the six Spitzer photometric channels, corresponding to the red and 
green curves, respectively. 



telescope (Borucki et al. 2009). However, modeling the 
visible flux in the Kepler bandpass introduces another 
five free parameters in terms of the dominant opacity 
sources in the visible - the atmospheric concentrations 
of TiO, VO, Na, K, and a prescription for scattering. 
These five parameters, which arc largely decoupled from 
the opacity sources in the IR, allow a large degree of flex- 
ibility in fitting the one Kepler observation. By exploring 
a preliminary range of values for the five parameters, we 
do find that some of our best-fit non-inversion models for 
the Spitzer data are also able to fit the Kepler point. In 
an earlier work, we found that the Kepler point could also 
be fit with models with thermal inversions (Christiansen 
et al, 2010). 

4.3. TrES-4 

TrES-4 is the second hottest planet in our sample, after 
HAT-P-7b, and is highly favored to host a thermal inver- 
sion, on theoretical grounds (Fortney et al. 2008). We 
use the five Spitzer photometric observations of Knutson 
et al. (2009), which have been previously reported as 
evidence for a thermal inversion in TrES-4. 

We find that existing observations of TrES-4 can be ex- 



plained almost equally well by models with and without 
thermal inversions, contrary to previous findings. Our 
fits for each case are better than £ 2 = 1, i.e within the 
l-cr error bars. The non-inversion P-T profiles fitting 
the data at the £ 2 < 1 and £ 2 < 2 levels are shown in 
Figure [31 

As shown in Figure [4j non-inversion models can fit the 
observations of TrES-4 to within £ 2 < 1 . for a wide range 
of CO and CH4 abundances, including the values at TE Q , 
represented by fco = 1 an d /ch 4 = 1- However, we find 
that the best-fitting models (within £ 2 < 1) without in- 
versions require a high concentration of CO2; > 10 -4 , for 
TE Q concentrations of the remaining molecules like H2O 
(e.g. Figure [5J. This requirement arises from the non- 
detection of planet flux observed in the 16 fim channel. 
Since the dominant contribution to this channel comes 
from the CO2 feature at 15 ^tm, a non-detection of flux 
indicates substantial absorption due to CO2 and hence 
the high CO2 requirement. At the high temperatures 
of TrES-4, a CO2 concentration of 10~ 4 is not feasible 
via equilibrium chemistry at solar abundances. However, 
CO2 concentrations as high as ~ 10~ 4 are, in principle, 
feasible for a high metallicity, over 30 x solar, (Zahnle ct 
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TABLE 3 

Chemical compositions and day-night redistribution of sample spectra 3. 





H a O 


CO 


CH 4 


co 2 


C/O 


V 


e 


HD 209458b 
I 

NI 


6 x lO" 06 (0.006) 
4 x 10~ 05 (0.04) 


3 x 10- 02 (50) 
8 x lO" 06 (0.01) 


4 x 10~ 08 (0.001) 
9 x lO" 06 (0.8) 


6 x lO" 09 (0.006) 
5 x lO" 11 (5 x lO" 05 ) 


1.00 
0.36 


0.38 
0.56 


1.60 
2.75 


HAT-P-7b 
I 

NI 


7 x 10-° 3 (7) 
2 x 10~ 05 (0.02) 


6 x 10-° 4 (0.8) 
5 x 10-° 4 (0.63) 


8 x lO" 07 (0.1) 
2 x lO" 11 (0.01) 


1 x 10~ 05 (10) 
8 x lO" 05 (80) 


0.07 
0.85 


0.59 
0.63 


0.89 
3.0 


TrES-4 
I 

NI 


8 x 10~ 03 (8) 
5 x lO" 03 (5) 


4 x 10~ 04 (0.6) 
4 x 10" 04 (0.6) 


5 x 10~ 08 (1) 
2 x 10" 9 (0.2) 


2 x 10~ 08 (0.02) 
2 x lO" 05 (23) 


0.05 
0.09 


0.40 
0.41 


1.76 
1.80 


TrES-2 
I 

NI 


2 x 10~ 03 (2) 
1 x 10~ 03 (1) 


1 x lO" 03 (2) 
4 x 10~ 04 (0.6) 


1 x 10~ 05 (0.3) 

2 x 10" 08 (0.7) 


1 x 10~ 08 (0.01) 

2 x 10~ 08 (0.02) 


0.44 
0.23 


0.72 
0.46 


0.54 
0.98 



a The two rows for each system correspond to the two model spectra for each system presented in Figure [6] "I" and "NI" 
correspond to the models with and without a thermal inversion, respectively, in Figure [6] For each molecule, the mixing ratio 
averaged over all the layers is reported, along with the deviation from TEq shown in parentheses. 



e < 2. 

The constraints on the atmospheric composition and 
day-night redistribution are shown in Tables Q] and [21 
and in Figures 2] and [5] . For each scenario, i.e., with 
or without a thermal inversion, there is a wide range in 
molecular composition that can explain the data, includ- 
ing that of thermochemical equilibrium with solar abun- 
dances (i.e TE©). The composition in each scenario is 
practically unconstrained even at the £ 2 < 1 level. It fol- 
lows that the C/O ratio is also unconstrained by the data, 
spanning a rather unphysical range of 5 x 10~ 5 - 4 x 10 3 
even at the £ 2 = 1 level, for the non-inversion models, 
for example. The large range in C/O is a consequence of 
allowing the molecular mixing ratios to vary arbitrarily. 
Nevertheless, it does indicate that best-fit solutions with 
no thermal inversions can be found at the £ 2 < 1 level 
with very plausible molecular concentrations (sec § 15. 2[) , 
and both carbon dominated and oxygen dominated at- 
mospheres are allowed by the observations. Similarly, 
the day-night redistribution is also unconstrained by the 
data, at the £ 2 = 1 level. 

Based on the existing Spitzer IRAC observations, 
therefore, our results show that a thermal inversion can- 
not be inferred in the day-side atmosphere of TrES-2. 
The weak constraints on the atmosphere of TrES-2 is 
evident from the data. The flux ratio in the 4.5 fj,m 
IRAC channel is noticeably higher than that in the 3.6 
/zm channel, hinting at a possible emission feature due 
to a thermal inversion. However, the flux ratio in the 5.8 
fj,m channel is noticeably lower than that in the 8 /im 
channel allowing for H2O absorption, and hence the lack 
of a thermal inversion. Thus, while a model with thermal 
inversion can explain the data, a thermal inversion is not 
required by existing IRAC observations of TrES-2. Two 
best-fit models for TrES-2, with and without a thermal 
inversion, are shown in Figure [6j and the corresponding 
compositions are shown in Table [3] The models for both 
the cases shown here span a rather plausible range of 
chemical compositions. 

5. SUMMARY AND DISCUSSION 

We have investigated the question of whether thermal 
inversions can be robustly inferred from existing Spitzer 



al. 2009; Madhusudhan & Seager, 2010). 

Non-inversion model fits at the £ 2 < 2 level, on the 
other hand, do allow for a very plausible set of chemical 
compositions. Figure [5] shows a model spectrum without 
a thermal inversion (in red) and having a chemical com- 
position close to TE (Table [3] shows the composition). 
Also shown for reference is a model with a thermal in- 
version. The constraints on the molecular abundances 
are shown in Table [1] and Table [2] for models without 
and with thermal inversions, respectively. We conclude 
that the observations place almost no constraints on the 
molecular concentrations in the atmosphere of TrES-4, in 
cither scenario. While the wide range in allowed chemi- 
cal compositions is seemingly implausible, it docs allow 
for a significant population of non-inversion models that 
are physically plausible (sec l5.2p . We find this evidence 
enough to conclude that there is no sure sign of a thermal 
inversion in TrES-4, given current data. 

4.4. TrES-2 

Our results show that the observations of TrES-2 can 
be explained to a high degree of fit by models both with 
and without thermal inversions. This general conclusion 
is similar to that of O'Donovan et al. (2010), where we 
first reported that the data could be fit by models with 
and without thermal inversions. In that study our non- 
inversion models fitting the data required a CO abun- 
dance that was lower than TE Q by about two orders of 
magnitude, whereas the best-fit inversion models allowed 
TEq composition. However, with the new parameter ex- 
ploration routine, in the present study we have been able 
to explore regions of parameter space well beyond what 
we could pursue in O'Donovan et al. (2010). While our 
present results for inversion models agree with our pre- 
vious findings, our results for non-inversion models go 
beyond our findings in O'Donovan et al. (2010). 

The best-fit solutions with no thermal inversions in the 
present study span a wide range in chemical composition, 
including that of TE Q . Our results show that the P-T 
profile is mostly unconstrained by the data, resulting in 
a large region of parameter space that can explain the 
observations even at the £ 2 = 1 level. Figure [3] shows the 
non-inversion P-T profiles corresponding to £ 2 < 1 and 
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photometric observations of thermal emission from hot 
Jupiter atmospheres. We addressed this objective by 
thoroughly exploring the parameter space of models with 
and without thermal inversions for a sample of four hot 
Jupiters. We considered four systems which have Spitzer 
observations at four or more wavelengths, and are highly 
irradiated, so that they are theoretically favored to host 
thermal inversions (Hubeny et al. 2003; Fortney et al. 
2008): HD 209458b, HAT-P-7b, TrES-4, and TrES-2. 
Furthermore, the observations considered have also been 
previously reported to be consistent with thermal inver- 
sions in the corresponding systems, albeit less robustly. 
In this work, we addressed to what level of statistical 
significance and physical plausibility thermal inversions 
can be inferred in each of these systems. 

Our primary finding is that a detailed exploration of 
the model parameter space is necessary to make robust 
inferences of thermal inversions in cxoplanctary atmo- 
spheres. We find that the observations of TrES-4 and 
TrES-2 can be explained by models with and without 
thermal inversions, and with physically plausible chem- 
ical compositions, at the £ 2 < 2 and £ 2 < 1 levels, 
respectively. This finding is in contrast to the find- 
ings of Knutson et al. (2009) and Spiegel & Burrows 
(2010) who, based on forward models of Burrows et al. 
(2008), have suggested the requirement of thermal in- 
version from the same data sets. The difference in con- 
clusions results from the different modeling schemes; the 
forward models of Burrows et al (2008) assume layer- 
by-layer radiative equilibrium and chemical equilibrium, 
and a parametrization for day-night energy redistribu- 
tion. On the other hand, for a given dataset, we ex- 
plore the space of temperature profile and composition, 
without the requirement of layer-by-layer radiative equi- 
librium or chemical equilibrium, but still imposing the 
strict constraint of global energy balance. We note that 
our best-fitting inversion solutions for TrES-2 and TrES- 
4 do include a wide range of profiles with thermal in- 
versions, which likely encompass the profiles of Knutson 
et al. (2008) and Spiegel & Burrows (2010). Another 
interesting result is that if TrES-4 indeed does not host 
a thermal inversion, the best-fitting non-inversion solu- 
tions require a large CO2 mixing ratio (> 10~ 5 ), which 
might be an indication of enhanced metallicity in TrES-4 
(Zahnle et al. 2009; Madhusudhan & Seager, 2010). 

For HD 209458b and HAT-P-7b, we find that the ob- 
servations cannot be explained without thermal inver- 
sions to within a £ 2 < 3, i.e to within 1.7a observational 
uncertainties, for any plausible composition. Any better 
fit would require substantial enhancements in methane 
and depletion of CO, which is implausible at the high 
temperatures in the systems considered. Our inference of 
thermal inversions in HD 209458b and HAT-P-7b is con- 
sistent with previous findings of Burrows et al. (2008), 
Knutson et al. (2008), Christiansen et al. (2010), and 
Spiegel & Burrows (2010). Our results show that a de- 
tailed exploration of the model parameter space and an 
accurate assessment of the observational errors is essen- 
tial to robustly infer thermal inversions based on existing 
photometric observations. For example, if one considers 
the 2-it error bars on the data, thermal inversions may 
not be required even for compositions in chemical equi- 
librium. 



5.1. Thermal Inversions or Not? 

Whether or not the observations considered in this 
work can be explained without a thermal inversion de- 
pends on what level of fit, and physical plausibility of 
models, one is willing to consider. If we do not consider 
the physical plausibility of the best-fit models, the ob- 
servations of all the four hot Jupiters can be explained 
without thermal inversions to within the 1.5 a error bars, 
i.e £ 2 < 2.25. On the other hand, if we consider only 
models fitting within the la error bars, and/or enforce 
arguments of physical plausibility (see ? !5.2l below). HAT- 
P-7b and HD 209458b emerge as likely to host thermal 
inversions. The inference of a thermal inversion is, there- 
fore, very sensitive to the reported observational uncer- 
tainties. 

The inference of a thermal inversion can also be sensi- 
tive to one data point over others. In HD 209458b, for 
instance, several of the six observations can be fit at the 
~ la error bars by non-inversion models (e.g., Figure [6]). 
A large contribution to the £ 2 comes predominantly from 
the 5.8 fim point which is fit only at greater than 2a. 
Thus, for these models it is only the 5.8 (im IRAC point 
which guides any inference wc make about thermal inver- 
sions. Similarly, for the model spectra shown for HAT- 
P-7b, the dominant contribution to £ 2 comes from the 
8 /zm IRAC point. Therefore, any inference of thermal 
inversions can be highly sensitive to the reported obser- 
vational uncertainties in a single channel, which varies 
on a case by case basis. 

At the level of current observations, our results poten- 
tially deviate from theory. Given that the atmospheres 
of TrES-4 and TrES-2 can be explained by models with- 
out thermal inversions, it is possible that these systems 
do not host thermal inversions. If that happens to be 
the case, the results are in contrast to theoretical predic- 
tions, as both TrES-4 and TrES-2 have higher levels of 
incident star flux as compared to HD 209458b (Fortney 
et al. 2008), and hence are more likely to host thermal 
inversions. Nevertheless, since the observations for these 
systems are consistent with models both with and with- 
out thermal inversions, the only conclusion is that it is 
too early to claim thermal inversions in these systems, 
contrary to some previous studies. 

5.2. Plausibility of Models 

An important point concerns the physical plausibil- 
ity of non-inversion models fitting the observations. As 
shown in Figure HI the best-fit non-inversion models for 
HD 209458b and HAT-P-7b require substantial enhance- 
ment of CH4 as compared to CO. However, as explained 
in § 12. 3[ in the very hot atmospheres of these planets, CO 
is expected to be the dominant molecule, based on atmo- 
spheric chemistry. We currently do not have a physically 
plausible explanation for such CH4 enhancement at the 
expense of CO in a very-hot atmosphere. 

Our best-fit models explore an unrestricted range of 
atmospheric compositions. In trying to conduct an unbi- 
ased exploration of the parameter space, we have allowed 
for all the molecules to vary over a large range of values, 
that might be seemingly unphysical. For example, in 
TrES-2 the £ 2 = 1 limits extend to mixing ratios as high 
as 0.1 for H2O, CO and CH4. Such high mixing ratios in- 
dicate extreme metallicities that are too high to be plau- 
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sible, although not impossible. For reference, solar abun- 
dances have number fractions of C and O at 3.3 x 10 -4 
and 7.8 x 10 -4 , respectively (Anders & Grevesse, 1989; 
used in equilibrium chemistry calculations of Burrows & 
Sharp, 1999; but c.f. Allendc Prieto et al., 2002; As- 
plund et al., 2005, for recent values which are lower by 
a factor of - 2). Similarly, for HAT-P-7b and TrES-2 
the C/O ratios span many orders of magnitude reaching 
values as high as 10 4 , which are also manifestly unphysi- 
cal. For instance, the solar C/O ratio is ~ 0.5. However, 
the best-fitting models, with and/or without inversions, 
for all the systems do allow C/O ratios in the plausible 
range of 0.1 - 1. 

Finally, we have not explored the realm of drastically 
inhomogeneous models - those models where the mixing 
ratio of a molecule could be very different in different lay- 
ers of the atmosphere. It is understandable that photo- 
chemistry and vertical mixing can deviate molecular mix- 
ing ratios away from equilibrium (e.g. Line et al. 2010; 
Madhusudhan & Seager, 2010). We expect that our pre- 
scription for molecular species, which is parametrized as 
deviations from chemical equilibrium, spans the space of 
possible deviations. However, we have not considered 
arbitrarily populating different layers of the atmosphere 
with different species. Adhoc filling of the layers with 
specific molecules might allow a large, albeit unphysical, 
degree of freedom in fitting the observations. 

5.3. Future Observations to Resolve the Degeneracy 

In this work, we have addressed the apparent degen- 
eracy between atmospheric composition and thermal in- 
versions in hot Jupiter atmospheres. The large number 
of model parameters allows the freedom to fit the lim- 
ited observations of some atmospheres in any scenario, i.e 
with or without inversions; although, at different levels of 
fit. Future developments in observations and theory are 
needed to break the apparent degeneracies. Theoretical 
efforts are needed to put limits of physical plausibility on 
the atmospheric composition and temperature structure. 
Such limits, for example, might exclude many of the non- 
inversion models that fit the observations considered in 
this work. 

New observations are important for better constraints 
on thermal inversions. In the near future, multiple obser- 
vations of thermal emission from transiting hot Jupitcrs 
in the near-IR, from ground and with HST, along with 
existing Spitzer data can help constrain models to a good 
extent, on a case-by-case basis. For instance, the two 
model spectra of HAT-P-7b in Figure|6]show only modest 
differences in the Spitzer bandpasses. However, the spec- 
tra are markedly different in the near-IR, especially in 
the continua between the molecular features, which can 
potentially be observed from ground, e.g., in the J, H . 
and K bands. Additional constraints can be also placed 
by observations within the molecular features, which are 
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possible with space-based observations (e.g. with HST), 
except for molecules like methane where ground-based 
observations might also be feasible. Near-IR observa- 
tions have been reported for a few systems to date, not 
particularly constraining thermal inversions (e.g. Swain 
et al. 2009; Croll et al. 2010). However, multi-band 
near-IR photometry and/or spectroscopy can prove to 
be a rich resource for targeted searches for thermal in- 
versions in exoplanetary atmospheres. Targets can be se- 
lected based on constraints from already existing Spitzer 
observations, irradiation levels, and sensitivity of a given 
instrument to the planet-star flux contrasts. The several 
near-IR bandpasses mentioned above are currently ripe 
for this purpose. In the long run, high resolution spectra 
with the James Webb Space Telescope will have the sen- 
sitivity to conclusively identify the presence of thermal 
inversions based on spectrally resolved emission features. 

An important point concerns the opportunity to ob- 
serve in the two IRAC channels on Warm Spitzer. While 
it is true that robust inferences of thermal inversions 
could not be made even with four Spitzer points in sev- 
eral known cases, a large difference between the 3.6 fiva. 
and 4.5 /zm channels can place stringent constraints on 
the existence of thermal inversions. A very large flux 
excess in the 4.5 /zm channel over the 3.6 /zm channel 
is highly indicative of a thermal inversion; although, it 
also depends on the irradiation level of the planet which 
governs the blackbody continuum. A large excess in the 
3.6 /im channel, on the other hand, is almost a sure sign 
of no thermal inversion, as in the cases of HD 189733b 
(Madhusudhan & Seager, 2009) and GJ 436b (Stevenson 
et al. 2010; Madhusudhan & Seager, 2010). Thus, obser- 
vations of hot Jupiters with Warm Spitzer would likely 
be able to identify the extreme cases of systems with or 
without thermal inversions. 

Tremendous progress has been made in the last decade 
in our understanding of exoplanetary atmospheres. At 
the same time, recent and current observations allow us 
a chance to revisit previous interpretations made with 
limited observations. It is now upon us to evaluate all 
the theoretical options and observational uncertainties so 
as to determine a framework in which to interpret obser- 
vations. Judicious target-selection and efficient planning 
of future observations, from ground and from space, will 
be critical to characterizing the atmospheres of the grow- 
ing number of transiting exoplanets. 
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